Temperature and Magnetic Field Dependencies 
in a Model of Diamagnetic Hysteresis in Beryllium 
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The model of diamagnetic hysteresis loop of strongly correlated electron gas at the conditions 
of dHvA effect is developed. It is shown, that in the framework of Rayleigh theory for magnetiza- 
tion loop, the coercive force and remnant magnetization in every period of dHvA oscillations are 
characterized by strong dependencies on temperature, magnetic field and Dingle temperature. 

PACS numbers: 75.20.En, 75.60.Ch, 71.10.Ca, 71.70.Di, 71.25.-s; 71.25. He; 75.40-s; 75.40.Cx. 



I. INTRODUCTION 



At high magnetic field and low temperature the 
strongly correlated electron gas is characterized by insta- 
bility due to magnetic interaction between electrons Q, 
which results in phase transition from homogeneous state 
into inhomogeneous one and stratification of the sample 
into diamagnetic domains, e.g. Condon domains (CDs), 
with alternating from one domain to another diamagnetic 
moment. This instability of an electron gas is called dia- 
magnetic phase transition and extensively studied, both 
theoretically and experimentally The diamagnetic 

domains were observed in Be by magnetization measure- 
ments 0; in Ag by nuclear magnetic resonance (NMR) 
and recently by Hall probe spectroscopy Q|. They 
have been observed in Be, Sn, Al and Pb by muon spin- 
rotation (juSR) spectroscopy @. Usually, for plate-like 
samples, this non-uniform phase is realized as periodic 
domain structure with alternating from one do- 

main to another diamagnetic moment. Magnetic inter- 
action of strongly correlated electrons gas results in non- 
linear relationship between local diamagnetic moments 
and external magnetic field. It was shown Q that in 
CD phase the local diamagnetic moments, which con- 
tribute to the process of magnetization, are characterized 
by strong temperature and magnetic field dependencies. 

The existence of metastable states for local diamag- 
netic moments in each period of dHvA oscillations can re- 
sults in hysteresis, e. g. irreversible changes in magneti- 
zation due to change in external magnetic field, similarly 
to magnetic materials of spin origin. The microscopic 
mechanism of magnetization loop, relates to defects of 
crystal lattice, intrinsic stress and elastic deformations. 
The impurities, dislocations and other defects of crystal 
lattice cause an internal friction which opposes to the 
wall displacement. The domain walls are pinned by such 
defects. Unpinning is realized by small dissipative jumps 
through the increasing of the value of applied magnetic 
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field. 

Although, the possibilities of irreversible behavior of 
magnetization curve in diamagnetic materials have been 
discussed earlier H, only recent experimental investiga- 
tion of CDs in Be Q reveals hysteresis loop in the mag- 
netization process. The hysteresis magnetization curve 
of Be @ was measured directly by Hall probes in dc field 
and standard ac method with different modulation lev- 
els, frequencies and magnetic field ramp rates, and the 
shape of magnetization curve was reconstructed by sev- 
eral higher harmonics of the ac pickup voltage. Detec- 
tion of the hysteresis in Be substantiates a far reach- 
ing analogy between an electron gas, which has unstable, 
metastable and stable states within one dHvA period, 
and the liquid-gasphase coexistence and transformation. 
Recent progress 0, @ m experiment on observation of 
Condon domain structure in Ag and Be, in particular, 
the detection of diamagnetic hysteresis, provides a nat- 
ural stimulus towards a more detailed understanding of 
the properties of strongly correlated electron gas in the 
conditions of dHvA effect. CDs in a 3d electron gas are 
the only type of magnetic structures, for which the mag- 
netization process has not been considered so far. 

In present paper we develop the model of diamagnetic 
hysteresis, which fits the recent data on observation of 
diamagnetic hysteresis in Be p|. The hysteresis loop ap- 
pears periodically in every period of dHvA oscillations, 
when [Iqx > 1 ( x lS differential magnetic susceptibility, 
M is the oscillatory part of the magnetization, and B 
is magnetic induction []]). There are three main con- 
tributions into the process of magnetization of normal 
metals in the condition of strong dHvA effect which have 
to be distinguished, the dependence of local diamagnetic 
moments on external magnetic field through the reduced 
amplitude of dHvA oscillations 0, the domain wall dis- 
placement and nucleation of the domains. The rotation 
processes, that are intrinsic property of magnetization 
loop in magnetic materials of spin origin, can be ne- 
glected in rearrangement of CD structure at least for 
metals with relatively simple Fermi surfaces pj due to 
one-dimensional character of the diamagnetic moments. 
We show that the remnant magnetization and coerciv- 
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ity, as well as the hysteresis loop square, which are the 
fundamental properties of the magnetization curve, have 
strong dependence on temperature, magnetic field and 
Dingle temperature. 

The paper is organized as follows. Section[I]is an intro- 
duction. In Section [H] we consider the model of diamag- 
netic hysteresis, based on Rayleigh model in the frame- 
work of Lifshitz-Kosevich-Shoenberg approximation. In 
Section ITTT1 we calculate the temperature and magnetic 
field dependence of remnant magnetization and coerciv- 
ity Section IV is conclusions. The last section contains 
acknowledgements . 



II. MODEL 

The oscillator part of the thermodynamic potential 
density in the Lifschitz-Kosevich-Shoenberg formalism 
can be written in reduced form with taking into account 
the shape sample effects: 



n = 



Airk 2 



a cos b + 2° 2 (1 ~ n ) sm2 ^ 



(1) 



where b = k{B — hqH) = k[fioh ex + (1 — n)M] = 
x + (l — n)y, [i H is the magnetic field inside the material 
at the center of dHvA period, k = 2ttF/ (^ H) 2 , F is the 
fundamental oscillation frequency), H ex is external mag- 
netic field, h ex = H ex — H is the small increment to the 
magnetic field H , n is the demagnetization factor. All the 
components of vectors are taken along the direction of the 
magnetic induction. The quality a = AirkA = nodM/dB 
is the reduced amplitude of the dHvA oscillations (the 
ratio between the amplitude of oscillations A and their 
period AH). 

In the first harmonic approximation the magnetization 
is found from the implicit equation of state 0: 



y 



a sm \x 



(1 - n)y], 



(2) 



The validity of thermodynamic potential density fl QJ, 
and the expression for the reduced magnetization y (J2J, 
derived from Eq. is restricted by applicability to ho- 
mogeneous phase only, where the existence of demagneti- 
zation coefficient n is justified. In the conditions of strong 
magnetic interaction between electrons, when 



small in comparing to the magneto-static energy of ho- 
mogeneous state. In the domain state, when the reduced 
amplitude of dHvA oscillations satisfies to the condition 
PJl. the demagnetization factor n, which is characteristic 
of the uniformly magnetized sample, is replaced by the 
coefficient a |7(. This coefficient depends on magnetic 
field, temperature, impurities and geometry of the sam- 
ple and takes into account the magneto-static energy for 
given type of domain structure. Typically for plate-like 
sample, a << 1, independently on the types of domain 
structures [To| and important mainly in calculation of 
the steady-state period 2D of the domain structure by 
minimization of the free energy density, containing two 
terms, magneto-static energy and the surface energy of 
the domain walls 0I"E3- For observed CDs structures 
0,0, the inequality D/L « 1, where L is the width 
of the sample, is fulfilled and a ~ D / L <<1. 
The equation 



a(p H,T,T D ) = 1, 



(4) 



defines the critical surface in three dimensions [IqH — T — 
Td ■ Above this surface the uniform diamagnetic phase 
exists, but below it, the CD phase can appear in the part 
of every period of the dHvA oscillations. 

For proper calculation of the reduced amplitude of the 
dHvA oscillations a and constructing the hysteresis loop, 
the correct topology of Fermi surface has to be taken into 
account. In case of Be the standard procedure of calcu- 
lation of the amplitude of dHvA oscillations 0, based 
on series expansions near the extreme cross-sections, re- 
quires considerable modification due to negligible small 
curvature of extreme cross-sections of the 'cigar'-like part 
of the Fermi surface, contributed to the dHvA oscilla- 
tions 01 • The Fermi surface of Be consists of the second 
zone monster ('coronet') and the third zone 'cigar'. It 
is well-established that the dHvA oscillations originates 
from the three maximum cross-sections of 'cigar' ('waist' 
and 'hips') which are characterized by a very small cur- 
vature. In particular, the small curvature of the cylinder 
like Fermi surface explains relatively high amplitude of 
dHvA oscillations in Be at magnetic field H applied par- 
allel to the hexagonal axis, e.g. H || [0001], in comparing 
to Ag 0,^3) where the small amplitude of dHvA oscil- 
lations is explained in the framework of spherical Fermi 
surface 0. In the framework of the model of slightly 
corrugated cylinder like Fermi surface of Be 0] , the re- 
duced amplitude of the dHvA oscillations takes a form 



a(fx H,T,T D )>l, (3) 

a state of lower thermodynamic potential is achieved over 
part of dHvA oscillation cycle by the formation of CD 
structure, for which the local value of magnetization al- 
ternates in sign from one domain to the next, resulting 
in essential reduce of magneto-static energy. The usu- 
ally observed diamagnetic domain structure is of stripe- 
domain type 0,0, and the contribution of magneto- 
static energy into the free energy density is negligibly 



a= 5 y' ' \QbioH)\R(noH,T,T D ), (5) 

where R(fioH 7 T, Tp) is a reduction factor 0, which 
takes into account the influence of the temperature T and 
impurities on the amplitude of the dHvA oscillations. 
The complex function 

Q = Qi+jQ2=\Q\exp (jV), 
Qi =3/?cosA + 2(1-2/3)J (A), Q 2 =/3sinA (6) 
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where A = 1 2 Aq£,, depends on applied magnetic field 
through parameter I 2 = ch/e^H and describes the beat- 
ings in the amplitude of dHvA oscillations resulting from 
two different extreme cross-section of Fermi surface of Be 
[l2| . Here, k = 0.217 A° _1 p| is the distance in recipro- 
cal space between two extreme cross sections of 'cigar', 
Aq is average cross section area, Jo (A) is Bessel function 
of the first order and [3 is the half-width of the cylinder 
like extreme cross sections. 

The Eqs. (@J,© define in explicit form the critical tem- 
perature for diamagnetic phase transition as a function 
of magnetic field fJ-oH and Dingle temperature Tp. It 
allows us to construct the phase diagrams for Be and 
investigate the properties of the diamagnetic hysteresis. 

III. RESULTS AND DISCUSSIONS 

At the condition a(noH, T, To) > 1 the sample is di- 
vided into domains with up and down magnetization. 
In the center of the dHvA period h = the domains 
are characterized by the same width. An external mag- 
netic field removes the equivalence of the states, the en- 
ergy balance is altered and rearrangements of the domain 
structure take place. For CD structure this process is 
realized not only through the motion of domain walls, 
which results in corresponding changes of domain vol- 
umes, and nucleation of the domains, but also through 
magnetic field dependence of the local magnetization in 
every adjacent domains [7J. To understand the magneti- 
zation process in CD phase it is convenient to define the 
approximate anhysteretic, or pinning-free, magnetization 

i anhys) = \{l-^)y- + \{l + (7) 

where the function y± — y±{x) = —y T (—x) is the so- 
lution of the equation y± = a sin {x + y±) at the range 
\x\ < A < k ■ AH/ 2 in every period AH of the dHvA 
oscillations and describes magnetic field dependence of 
the local diamagnetic moments of the up (+) and down 
(— ) domains (Fig. [IJ. In Eq. (JJJ we assumed that the 
volume fractions of the up and down domains are linear 
functions of the increment of the reduced magnetic field 
in every dHvA period. This assumption was justified by 
NMR experiments in Ag and /iSR spectroscopy in Be 
HI- Both experiments reveal the domain structure by 
the appearance of a doublet corresponding to two domain 
magnetic inductions with line intensities, that show the 
linear dependence on the external magnetic field of the 
volume fractions occupied by the neighboring domains. 
The range of existence of the CDs 2A is defined by |7J 

A = ^Ja 2 - 1 - cos" 1 - > 0, a>l. (8) 
a 

It follows from Eq. J7J), that the anhysteretic magnetiza- 
tion in every period of the dHvA oscillations is a function 
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FIG. 1: Diamagnetic magnetization loop y — y(x; a) versus 
magnetic field increment x is plotted in one period of dHvA 
oscillations. The local diamagnetic moments y± = y±(x;a) 
reproduced from contribute to the hysteresis curve y a ,d = 
y a ,d{x;a). Three fundamental parameters of the hysteresis 
loop, e.g. coercivity x c , remnant magnetization y r and satu- 
ration magnetization y s are shown. The dot lines correspond 
to the range of existence of the hysteresis (2A). The exis- 
tence of the hysteresis is restricted by 2A < 2ir. The dash 
lines show the states which are not realized. In numerical 
calculations the value of a — 3 and n — 0.52 are used. 



of magnetic field H , temperature T and Dingle tempera- 
ture Tu- The expression for y( anh v s '> Eq. J7J) can be sim- 
plified for relatively high reduced amplitude of the dHvA 
oscillations, e.g. a>1. In this case y + « yo — (x/A)y , 
y- ~ y+ — 2?/o, where yo is the half-splitting in the cen- 
ter of the period, defined explicitly form by the equation 
yo / sm J/o = a HI, and the Eq. Q is replaced by the fol- 
lowing one 

y (anhys) „ ( 1 _ 1 ) |x| < A < 7T. (9) 

A 7T 

Although, the local diamagnetic moments y± = y±{x ex ) 
are characterized by average negative slope as the func- 
tions of magnetic field increment x ex 0, it follows from 
Eq. (|5J), that the anhysteretic curve has a positive slope 
in every period of the dHvA oscillations (Fig.QJ. We will 
see that the diamagnetic hysteresis is also characterized 
by the positive slope: in the process of the diamagnetic 
magnetization the volume fraction of the up-domains in- 
creases on expense of the volume fraction of the down- 
domains with increase of magnetic field increment x ex 



The real diamagnetic magnetization process in CD 
phase is characterized by hysteresis, or magnetization 
loop, which appears periodically in applied magnetic field 
in every period of dHvA oscillations, when the condition 
© is fulfilled. According to experimental data [4] the 
diamagnetic hysteresis is usually small, related to the do- 
main wall displacement and domain nucleation. Thus, it 
can be described in the framework of the Lord Rayleigh 
model with two parameters: initial susceptibility x an d 
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Rayleigh constant r] (see, e.g. [Tl|): 

y a ,d = ±y s + x{xT A) ± -77(0; T A) 2 , (10) 

where the upper (lower) sign corresponds to the ascend- 
ing (descending) magnetization in increasing (decreas- 
ing) magnetic field (Fig. The three fundamental 
parameters of the hysteresis loop: the remnant mag- 
netization y r = y r (a), the saturation magnetization y s 
and the coercivity x c = x a , which is the reverse field 
to bring the magnetization to zero from initial satura- 
tion, are characterized by strong dependence on the ap- 
plied magnetic field, temperature and Dingle tempera- 
ture through the reduced amplitude of the dHvA oscilla- 
tions a — a(fioH, T, Td). The saturation magnetization 
ys = ya can be calculated numerically from the explicit 
equation 



y s 



a sin [A + (1 — n)y s \. 



(11) 



The Eq. (|10() gives the following expressions for the co- 
ercivity factor x c and remnant magnetization y s : 



SA 2 



1-SA+ ^{1-SA) 2 + {SA) 2 
Vr = ^A 2 , 



SA 2 



(12) 
(13) 



where 6 = rj/x- In calculation of Eqs. i|12|) . (|13(l we used 
SA <^ I, which is justified by the observed narrow hys- 
teresis loop in Be |8|. In this limit, coercivity factor x c 
and remnant magnetization y ri being quadratic functions 
of the range of CD existence A ©, have the same tem- 
perature and magnetic field dependencies. The initial 
susceptibility \ can be calculated from measured rem- 
nant magnetization y r and coercivity factor x c according 
to the equation x = 2y r /x c . The hysteresis loop square 
which characterizes the losses in the process of diamag- 
netic magnetization, is defined as follows 



S= -77A 3 . 
3 ' 



(14) 



and also temperature and magnetic-field dependent. The 
loop square S Eq. (|14l) reaches its maximum S max — 
(4/3)7r 3 77 at a 3> 1, inasmuch as the range of the existence 
of the CDs 2A in every period of the dHvA oscillations is 
restricted by the period of the oscillations ( 2ir in relative 
units). Analysis of the Eqs. (|12|) - (|14fl with taking into 
consideration Eq. (jHJ shows that the parameters of the 
hysteresis loop x c , y r , S — ► at the point of diamagnetic 
phase transition, when a — ► 1. 

The magnetic field dependence of the coercivity x c — 
x c ([1qH) HI 21 1 at the Dingle temperature To = 2K of the 
experiment [8| and three different temperatures T from 
0.5K to 1.5K is represented in Fig. J5J. To construct 
the plots we used S = 0.1, which is justified by com- 
parison with the data |l3j . The function x c = x c (ijlqH) 
is plotted at the whole range of applied magnetic field 




3,5 4,0 4,5 
M H(T) 



FIG. 2: Magnetic field dependence of the reduced coercivity 
(fM)H) as a function of applied magnetic field at T = 
0.5K and Td = IK in the whole interval of CD existence in 
Be (a) and at three different temperatures T = 0.57^ (solid 
line), T — IK (broken line) and T = 1.5A" (dot line) (b) in 
the part of the magnetic field interval. For calculations the 
value of 5 = 0.1 is used. 



H from IT till 6T (Fig. where the non- uniform 

diamagnetic phase exists at given Dingle temperature 
Td = 2K The function x c — x c ((j,qH) has asymmet- 
ric multi-bell like shape, ty pica l for the phase diagrams 
of Be n H -T - T D [jj, [H. The maximums of the 
function x c = x c (hqH) appear periodically on the scale 
of inverse magnetic field with the period inversely pro- 
portional to the discrepancy AF — — F w of the two 
fundamental frequencies, F% and F w , corresponding to 
two extreme cross sections of cigar-like Fermi surface of 
Be: 



1 



1 

AF' 



(15) 



For F h = 970.9T and F w = 942. 2T we obtain 

AF = 28. IT . The increase of the temperature T re- 
sults in decrease of the maximums of the function x c = 
x c (^,qH) till its disappearing, when the point of phase 
transition is reached (a = 1). The position of the maxi- 
mums of the function x c — x c (fj,oH) does not depend on 
temperature in accordance with the phase diagrams for 
Be [l2| • The periodicity of the coercivity on reciprocal 
magnetic field is illustrated in Fig. Eta), where the co- 
ercive force X c = x c /k = A c (/ioF) measured in Gauss 
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is plotted as a function of (fioH)^ 1 at T = 1.3K and 
To = 1.8K. It follows from Fig.^a), that the preferable 
interval of applied magnetic field for measuring the coer- 
cive force is 2 — 4T , where the function X c can reach max- 
imum values from 1G to 2G at given values of T — 1.3K 
and T D = 1.8K. 

The increase of the Dingle temperature Tu — h/2nkBT 
which is inversely proportional to the scattering lifetime 
t of the conduction electrons, leads to the reduction of 
the amplitude of the dHvA oscillations, resulting in the 
phase transition into uniform diamagnetic phase. Thus, 
the range of the existence of the CD phase A Eq. (JSJ 
decreases (A — » 0, when a — > 1). Assuming that this 
contribution to the coercive force x c Eq. (|12Jl is dominant 
and neglecting an unessential change of the parameter 
5 = x/ 7 ?! one can come to the bells-like magnetic field 
dependencies of illustrated in Fig. [21 

The temperature dependencies of the reduced coer- 
civity factor x c — x c (T) and measured coercive force 
X c = X C (T) at /io-ff = 3.6T and different Dingle tem- 
peratures are shown in Fig. The circles in Fig. EJb) 
represent the results of measurement of the coercive force 
in Be 13]. The theoretical curve X c — X C {T) plot- 
ted a,tfi H = 3.6T and T D ps 1.8K fits the experimen- 
tal points. The reported Dingle temperature in ^| is 



2,0 
— 1,5 
X°1,0 
0,5 
0,0 
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A 
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T= 1.3 K 
T D =1.8K 




0,2 0,3 0,4 0,5 0,6 




FIG. 3: (a) Coercive force X c — x c /k — X c (/j,oH) as a func- 
tion of inverse magnetic field at T = 1.3K and Td = 1.8A", 
and (b) reduced coercivity 

Xc — Xc Oo-ff) at T — 1.3K and 
three different Dingle temperatures: To = 1-8 (solid line), 
T D = 2K (broken line) and T D = 2.2K (dot line). 



To = 2K. There is a good agreement between predicted 
temperature dependence of the coercive force and data 

m 

The observation of the hysteresis loop in Be |j] pro- 
vides the possibility of the verification of the developed 
theory of the diamagnetic hysteresis and evaluate pa- 
rameters of the loop. The periodical irreversible behav- 
ior of the magnetization process in the rod-like sample 
8 x 2 x 1mm 3 of Be is demonstrated at T — 1.3K in the 
whole range of applied magnetic field /io-ff ~ 2 — 6T, 
where the CD phase exists. The Dingle temperature 
is Td = IK. The diamagnetic phase diagrams esti- 
mated from Eq. 10} with taking into account the Eq. (J5J 
for the reduced amplitude of oscillations a, obtained in 
the framework of the slightly corrugated Fermi surface 
of Be 12J, are consistent with the data [T^- At 
[i^H = 3.6T the calculated period of the dHvA oscil- 
lations AH = 2n/k = (li H) 2 /F w 138G (F = 955T is 
average area of the extreme cross-sections of the Fermi 
surface sheet, relevant to dHvA oscillations Q) is very 
close to the observed period of the dHvA oscillations 
AH w 134 - 138G ||. According to the data @, the 
coercive force X c « 1.2G, the remnant magnetization 
Y r = y r /k 1.2G, the saturation magnetization Y s 
and the range of existence of the non-uniform phase 
2X\ can be evaluated as Y s w 30G and 2X& ~ 34G 
consequently. Thus, the experimental value of the re- 
duced parameter A^ exp ^ = X&k w 0.8 can be compared 
with the corresponding theoretical value of A Eq. © 
in the following way. At T = l.3K,T D = 2K and 
fioH = 3.6T, the reduced amplitude of the dHvA oscilla- 
tions according to the Eq. JSJ is a « 2, giving the value 
of A k 0.7 in accordance with the experimental value 
A (ex P ) = w o g ^ According to this calculation, 

at the conditions of the experiment , CD phase occu- 
pies w 0.25 part of the period. The data allow us to 
evaluate the parameters of the Rayleigh hysteresis loop, 
X and rj. First, we calculate y r = kY r w 0.053, that gives 
r] = 2y r /A 2 w 0.18, x = 2y r /x c w 2 and S = 0.09 « 0.1. 
Then, we can compare the experimental and theoretical 
values for saturation magnetization. According to the 
data yi exp) = kY s w 1.4. The theoretical value of 
the saturation magnetization y s can be calculated from 
Eq. (| 1 1 p . Using the value of the demagnetizing coefficient 
n w 0.05, calculated for prolate ellipsoid [l4|, a ~ 2 and 
A « 0.7 one can obtain y s w 1.6 in accordance with the 
experimental value of y e s xp = 1.4. 



IV. CONCLUSIONS 

Magnetic interaction of strongly correlated electron 
gas results in non-linear relationship between magneti- 
zation and magnetic field, which can give rise an effect 
of hysteresis in magnetization curve due to existence of 
metastable states for local diamagnetic moments in each 
period of dHvA oscillations. There are three main con- 
tributions into the process of magnetization of normal 
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0,15 




T(K) 

FIG. 4: (a) Temperature dependence of the reduced coercivity 
factor x c = Xc(T) and (b) coercive force X c = x c /k = X C (T) 
at /xo-ff = 3.6T and six different Dingle temperatures Tn : 
1 - 1.8if,2 - 1.9if,3 - 2K,4 - 2.LR r ,5 - 2.27^,6 - 2.3K. 
Close circles are referred from R.Kramer et al fl3|l . 

metals in the condition of strong dHvA effect: the do- 
main wall displacement, the nucleation of domains and 
the dependence of local diamagnetic moments on exter- 
nal magnetic field Q- 

We developed the model of diamagnetic hysteresis in 
CD phase in the framework of Rayleigh hysteresis loop. 
Rayleigh law for the magnetization process of spin origin 
is generally believed to be caused only by the wall dis- 
placement and in relatively small magnetic field. There 
are two circumstances that can justify the applicability of 
the Rayleigh model for magnetization process of strongly 
correlated electron gas. The first, although for observa- 
tion of the dHvA effect the high magnetic field in the scale 
of 1 — 10T is need, the interplay between uniform and non- 
uniform phases realized at the fine magnetic field scale 
of the magnetic field increment \h\ < AH = (^oH) 2 /F, 
e.g. ~ mT, in every period of the dHvA oscillations. 
And, the second, the magnetization process for electron 
gas in the conditions of dHvA effect is dominated by do- 
main wall motion which is similar to the magnetization 
process of soft thin magnetic films in the physics of mag- 



netism of spin origin. There is a good agreement be- 
tween the theoretical results obtained in the framework 
of the developed model of diamagnetic hysteresis and re- 
cent experimental results on observation of the hysteresis 
in Be |8(- The diamagnetic hysteresis is usually small, 
appears periodically in every period of dHvA oscillations 
with changing the value of applied magnetic field, and is 
described in the framework of the Rayleigh model with 
two fundamental parameters: initial susceptibility \ and 
Rayleigh constant r\. We calculate the parameters of the 
diamagnetic hysteresis from the data Q- The zero field 
saturation (remnant) magnetization y c , coercivity x c and 
hysteresis loop square S are characterized by strong de- 
pendence on temperature, magnetic field and impurity of 
the sample through the reduced amplitude of the dHvA 
oscillations a(fioH,T,TE>). The irreversible changes in 
magnetization are the result of existence of metastable 
states due to defects of crystal lattice, intrinsic stress 
and deformation. Generally, the diamagnetic hysteresis 
occurs because the domain boundaries do not return to 
their original positions when the external magnetic field 
approaches to its value corresponding to the centre of 
period of dHvA oscillations. At the initial permeabil- 
ity range the domain walls are displaced reversibly from 
their stable positions which essentially determined by the 
homogeneity of materials. Beyond the initial permeabil- 
ity range, the intensity of magnetization increases and 
the process is irreversible. The irreversible magnetization 
range is attained by irreversible displacement of domain 
walls. Although the rotation processes, that are intrin- 
sic property of magnetization loop in magnetic materials 
of spin origin, can be neglected in rearrangement of CD 
structure, the microscopic theory of diamagnetic loop has 
to include into consideration the possible bending of the 
domain wall due to magnetization current densities lo- 
calized in the domain wall close to the sample surface 
|9|. The irreversibility of the magnetization process is 
believed to be the intrinsic property of the CD structure 
for all normal metals. The present model of diamagnetic 
hysteresis can be applied for investigation of temperature 
and magnetic field dependencies of diamagnetic magne- 
tization not only in Be, but in other systems, where the 
CDs were observed. We hope that our theoretical results 
will motivate the experimentalists to carry out further 
experiments in investigation such an exotic phenomenon 
as diamagnetic hysteresis. 
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